Low-temperature nucleation in a kinetic Ising model with soft stochastic dynamics 
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We study low-temperature nucleation in kinetic Ising models by analytical and simulational meth- 
ods, confirming the general result for the average metastable lifetime, (r) — Aexp{/3T) (/3 — l/ksT) 
[E. Jordao Neves and R. H. Schonmann, Commun. Math. Phys. 137, 209 (1991)]. Contrary to com- 
mon belief, we find that both A and F depend significantly on the stochastic dynamic. In particular, 
for a "soft" dynamic, in which the effects of the interactions and the applied field factorize in the 
transition rates, F does not simply equal the energy barrier against nucleation, as it does for the 
standard Glauber dynamic, which does not have this factorization property. 

PACS numbers: 64.60.Qb 64.60.My 02.50.Ga 75.60.Jk 



Nucleation is fundamentally important in disciplines 
ranging from biochemistry [ij to earth sciences as- 
trophysics 0, and cosmology 0], and it has been stud- 
ied by kinetic Monte Carlo (MC) simulations in electro- 
chemistry ||5j|, materials science magnetism 0, and 
atmospheric science S], to mention just a few. However, 
many questions in nucleation theory are still unresolved, 
and recently there has been much interest in kinetic Ising 
systems as models for nucleation. In particular, much 
work has been done on their dynamical behavior at very 
low temperatures 0, E3, El El El El , where it is in- 
fluenced by lattice discreteness. It is then possible to 
calculate exactly both the shape of the critical nucleus 
(the saddle-point configuration) and the most probable 
path during a nucleation event. In a typical numerical 
experiment, the system is prepared in a metastable state 
with all spins positive in a negative applied field. During 
each MC step (MCS), a randomly chosen spin is fiipped 
with a configuration-dependent transition rate W that 
satisfies detailed balance, so that it will drive the sys- 
tem to thermodynamic equilibrium. The metastable life- 
time is measured as the average number of MCS until 
the magnetization reaches zero. In the regime of single- 
droplet decay studied here, the lifetime measured in MCS 
is independent of the system size H . In the low- 

temperature limit the lifetime has been rigorously shown 
to be li 



(r) = Ae/r. 



(1) 



Here the only dependence on the temperature T is 
through /3 = l/fceT, where fee is Boltzmann's constant 
(hereafter set equal to one). It is often assumed that P 
equals the energy difference between the saddle point and 
the metastable state ^1 Q| , independent of the specific 
stochastic dynamic. In this Letter we show that this is 
not always so. In particular, we describe two dynamics 



that both obey detailed balance but have different val- 
ues of P and A for all values of the applied field, despite 
having the same saddle-point configuration. 

At sufficiently low T, the saddle point was shown in 
Ref. |9J to be an ^ X (^ — 1) rectangle of overturned spins 
with a "knob" of one overturned spin on one of its long 
sides. The critical length I = [2J/\H\\ + 1 for all \H\ G 
(0,4), where \_x\ is the integer part of x. Here, J > 
is the nearest-neighbor interaction constant of the Ising 
model, which will henceforth be set to unity. The critical 
length thus changes discontinuously at values of such 
that 2/|iJ| is an integer. 

The square- lattice S = 1/2 Ising ferromagnet with 
unit interaction is defined by the Hamiltonian TL — 
-J2{a,i3)^acrp - Hj2a^a^ whcrc the Ising spins CTq = 
±1, H is the applied field, J2{a p) over all nearest- 
neighbor bonds on a square lattice, and 'Ylia runs over 
all lattice sites. When this system evolves under a 
continuous- time Glauber dynamic with spin-flip rate E3 



I^G = [1 + exp {fil^E)] 



(2) 



where /S.E is the energy change that would result from 
the flip, P in Eq. ^ is given by 
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(3) 



and from Ref. [ij A = ^Hard = 3/[8(^ - 1)] for aU \H\ < 
2. (See explanation of the subscript "Hard" below.) The 
interpretation of Pnard is indeed the energy difference 
between the saddle point and the metastable state. 

A characteristic feature of the Glauber dynamic is that 
it does not factorize into one part that depends only on 
the change in interaction energy, AEj, and another that 
depends only on the change in the field energy, AEh- 
Such transition rates are known as "hard" Il7|. 



TABLE I: Rates of flipping a spin in class m, pm, in the limit /3 oo. Here a = + ( — ) corresponds to a spin in the metastable 
(stable) direction, and A''+ is the number of its nearest-neighbor spins in the metastable direction. The analytic form of pm for 
the soft dynamic does not change with \H\. 
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Dynamics that do factorize this way are called "soft." An 
example is the soft Glauber dynamic [l^ . 



T^SG = [1 + exp iPAEj)]-^ [1 + exp {fSAEn)]- 



(4) 



In studies of field-driven Ising and solid-on-solid inter- 
faces ^3 it was recently shown that soft dynamics 
yield significantly different microscopic interface struc- 
tures and mobilities than hard dynamics. Here we show 
that also the low-temperature nucleation properties with 
the soft Glauber dynamic differ significantly from those 
with the hard Glauber dynamic. In particular, F is not 
simply the energy difference between the saddle point and 
the metastable state, and the prefactor A is also different. 

We obtain our results in three different ways. First, 
we calculate analytically by hand the first-passage time 
from the metastable state to an absorbing state just be- 
yond the saddle point in an approximation that the path 
in configuration space corresponds to a simple one-step 
Markov process |2flj . Second, we perform computer-aided 
analytical calculations u sing the technique of absorbing 
Markov chains (AMC) [ifl HJ, allowing for multiple 
branching paths and "blind alleys." Third, we perform 
simulations using the MC with AMC (MCAMC) tech- 
nique [lol E3 |. The first method provides the clearest 
physical insight, and for noninteger values of 2/\H\ the 
results are fully confirmed by the other two. 

The one-step Markov chain for 1 < |iJ| < 2 = 2) 
corresponds to the configurations labeled i — 0, ...,4 in 
Fig.n The label i gives the number of overturned spins, 
such that the starting configuration has i — 0, and the 
saddle point has i = i* = 3. In general the absorbing 
state is labeled / > i* > 1. The mean time spent in state 
i is hi. The rate at which the cluster grows from « to i -I- 1 
overturned spins is gi , and the rate with which it shrinks 
from i to « — 1 is Si. (Multiple spin fiips are negligibly 
rare in the zero-temperature limit These quantities 
satisfy the relation |20lEll23| 



(5) 



with boundary conditions S/ = Sq = 0. The number 
N of sites in the system represents the total probability 
current through the Markov chain . From Eq. lO we 
obtain hi recursively as = N/gi-i and 
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E — n 



(6) 



for < i < / — 2. Assuming that I > i* and the growth 
time for the supercritical droplet is negligible compared 
with the nucleation time fl^ , the lifetime (r) is 
the mean first-passage time to /, (t) = (t/) = J2i=o 
Grouping the terms according to the "unpaired" factors 
N/ gi+k in Eq. 10) then yields 
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(7) 



This result is general for any one-step Markov chain with 
absorption at /, regardless of the values of gi and Si 
|2(ij |. However, for the Ising model the transition rates 
are related by detailed balance as {si/nf)/{gi-i/nf_-^) — 
^l3(Ei-Ei-i) ^ where Ei is the energy of state i. The de- 
generacy factors and n^_j^ are the numbers of lattice 
sites at which a single spin fiip can shrink the cluster 
from « to i — 1 and analogously for growth from i — 1 to 
i, respectively. As a result, Eq. becomes 
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(8) 

In the limit f3 oo Eq. (jS)) is dominated by the term 
or terms with the largest exponential factor. Their selec- 
tion, which determines F and A, is described below, after 
we next find the spin- fiip rates in the different dynamics. 

For the square-lattice Ising system, the spins fall into 
10 classes |22|, determined by the spin value a {+ for the 
metastable direction and — for the stable direction) and 



2 



E^O E^=12-4\H\ E^=16-8\H\ 

90= ^Pl 91=4P2 02=4^2 93=^3 
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FIG. 1: The states in the one-step Markov chain of clusters of 
i = 0, 4 overturned spins, used to calculate the metastable 
lifetime (r) analytically by hand. The right-pointing arrows 
give the growth rates gi-i, and the left-pointing arrows give 
the shrinkage probabilities Si for i = 1, 2, and 3. The energies 
Ei (relative to the metastable state, i = 0) are given at the 
top of the figure for even i and at the bottom for odd i. 



= A + B + C + V. (10) 

Using p„i from Table^we identify the dominant terms in 
(T4), and from these we obtain A and F for both the soft 
and hard Glauber dynamics for all \H\ > 1. (Analogous 
calculations can be carried out for arbitrarily small \H\.) 

Soft dynamic: For 1 < \H\ < 2 the sum is dominated 
by the last term in C, yielding Asoft = 1/8 and Fgoft = 
16 - 4:\H\. For \H\ = 2 it is dominated by A and the last 
terms in B and C, yielding Agoft = 11/8 and Fsoft — 8. 
For \H\ > 2 it is dominated by A, yielding ^soft = 1 and 
Fsoft = 8. See Fig. El 

Hard dynamic: For 1 < \H\ < 2 the sum is dominated 
by the last terms in C and V, yielding Ajjard = 3/8 and 
Tnard = 16 — 6\H\. For \H\ = 2 it is dominated by A and 



the number of its nearest neighbors that point in the 
metastable direction. The low-temperature limits of the 
rates Pm for flipping a spin in class m [Eqs. ^ and Q], 
are shown in Tabled 

Figure shows a one-step Markov chain with I — A. 
For 1 < |iJ| < 2 the saddle-point configuration {i — i*) 
is the L-shapcd cluster with ^ = 2 (i* = 3). Among 
the dominant terms in Eq. ||SJ) is always the term with 
k = I = i* -1. For < |iJ| < 2, growth from i* - 1 to i* 
always involves adding a "knob" to one of the long sides 
of an ^ X (^ — 1) rectangle, such that f/i»_i = 2£p2. From 
TableHlwe see that = e-^^lffl^Hard p^^. 2 < < 

4, the saddle point is a single overturned spin {£ — 1), 
so i* = 1 and gi*-i — go — Npi. Again, the difference 
between F for the two dynamics is determined by the fact 
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This yields our main result: 
2|iJ| for0<|i7|<4. (9) 



We emphasize that Eq. © is valid as T ^ for all 
\H\ € (0,4), although the low-temperature regime will 
only be be reached at exceedingly low T as \H\ decreases. 
For > 4 the lifetime is the first-passage time to one 
overturned spin, so that (r) = (ti) = 1/pi, which yields 
Tsoft = 8 and Fnard — 0. Thus, in contrast to the hard 
dynamic, nucleation with the soft dynamic is always ac- 
tivated, even for infinitely strong fields. 

To obtain the prefactors Asoft and ^Hard in the two 
dynamics, we explicitly write out the four terms obtained 
from Eq. (jSJ for / = 4 with Si, gi-i, and Ei from Fig. ^ 



(^4) = - + -r- 



N - 



1 1 

Pi 4p2 
1 

+ — ( N 

P3 

N 
+ 4^^ 



N + e 
N 



/3(8-2|H|) 



„/3(12-4|ff|) 



2\H\) 



4|ff|) I _p/3(16-6|H|) 
^ 4 



20 r 



15- 








1 1 1 


1 






+ 


IHI = 1.9 








A 


IHI =2.0 










IHI =2.5 










Analytical 






V <7 ^ V ^ 







0.02 0.04 0.06 0.08 
T 



5- 



Soft, analytical 
Soft, simulated 
Hard, analytical 
Hard, simulated 



IHI 



1.5 



0.5 



(b) 



Soft, analytictil 
Soft, 1-step analytical 
Soft, 12/9 analytical 
Soft, simulated 
Hard, analytical 
Hard, 1-step analytical 
Htird, 12/9 analytical 
Hard, simulated 



IHI 



FIG. 2: Analytical and simulated results for the soft and 
hard Glauber dynamics for (a) F and (b) A. In the legends, 
"1-step analytical" refers to the one-step Markov-chain ap- 
proximation, "12/9 analytical" to the computer-aided AMC 
calculations with 12 transient and 9 absorbing configurations, 
and "analytical" to results that are identical for all the ana- 
lytical calculations. The results only differ for \H\ = 2. The 
inset in (a) shows the analytical (lines) and MC (data points) 
results for rin{r) vs T, from which F and A are obtained. 
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the last terms in B, C, and P, yielding Anard = 2 and 
Fnard = 4. For 2 < \H\ < 4 it is dominated by A and 
the last term in B, yielding Auard = 5/4 and Ffjard 

8 — 2\H\. (These results agree with corresponding ones 

in Refs. |l3Lll4l |.') For \H\ > 4 the system is unstable and 
(^) = (n) = A, yielding Fnard = and ^Hard = 2 for 
\H\ = 4 and ^Hard = 1 for \H\ > 4. See Fig.|21 

We further performed computer-aided analytic calcu- 
lations of (r) with Mathematica 25] by the AMC method 
[lol |2]| , using three different classifications of the config- 
urations: 12 transient and 9 absorbing states (denoted 
12/9), as well as 7/13 and 13/13. For noninteger 2/\H\ 
the results were identical to the 1-step approximation. 
However, for \H\ = 2, A (but not F) was found to depend 
slightly on the numbers of states included in the calcula- 
tion for both dynamics. With the numbers of states used, 
these differences were less than 0.5%. Specifically, 12/9 
yielded ^Hard = 78244/45597 w 1.7160 (1.764 by a differ- 
ent method in Ref. ^) and Agoft = 943/704 « 1.3395. 
See Fig. El 

Both sets of analytic results were checked by MC sim- 
ulations for both dynamics, using the MCAMC method 
[ifil E^ . The system size was L = 24, and 2000 es- 
capes were used (6000 for \H\ > 4). The parameters 
F and A were determined from weighted two-parameter 
linear least-squares fits to plots of Tln(r) vs T [inset in 
Fig-EIa)]- As seen in Fig. El the simulation results agree 
with the analytical results to within two standard errors. 

In conclusion, we have confirmed Eq. .9^ for the low- 
temperature metastable lifetime of a kinetic Ising model, 
using both analytical methods and MC simulations, find- 
ing both F and A to depend on the specific stochastic 
dynamic for all values of the applied field. For a soft 
Glauber dynamic [Eq. @], F does not equal the energy 
difference between the critical cluster and the metastable 
state for any value of the field and it also does not vanish 
in the strong- field limit, as it does for the conventional, 
hard Glauber dynamic [Eq. lEJ] . Thus, nucleation under 
the soft dynamic remains an activated process for ar- 
bitrarily strong fields. These results are consistent with 
recent studies of the microstructure and mobility of field- 
driven Ising and Solid-on-Solid interfaces [l^, UJ| • They 
indicate that great caution must be shown in formulating 
and interpreting stochastic models of physical systems, 
as even seemingly minor modifications of the transition 
probabilities can significantly affect the nucleation rates. 
It might thus be interesting to investigate the influence of 
the specific stochastic dynamic on dynamic phase transi- 
tions in kinetic Ising models [2^. We also note that, al- 
though our results are derived for a specific model system, 
qualitatively similar results should apply to kinetic MG 
simulations for nucleation in a wide range of scientific dis- 
ciplines. On the positive side, experimental observation 
of the field and temperature dependences of nucleation 
and growth could help devise correct stochastic models 
of nonequilibrium phenomena. 
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